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ABSTRACT 


Gihemevere active cevyices are included in an electronic circuit that 


is to be analyzed by a cont ee. appropriate models for these devices 
must be developed. A lumped large-signal dynamic model of the field- 
effect transistor (FET) is presented and the procedure for pointwise 
linearization of this model is described. This linearized model is 
suitable for use with the TRAC (Transient Radiation Analysis by Com- 
puter program) network analysis program. Implementation of this 
model using TRAC coding was demonstrated by programming an example 
circuit for each of two basic types of field-effect transistor. The 
performance of the model in simulating a basic pulse invertor circuit 
was compared with actual device behavior. Suggestions for extension 


of this work are included. 





TABLE OF CONTENTS 


I. INTRODUCTION SEI NN IIS a aa alata 
ILite BACKGROUND ------ nnn nen en en ee en en nr nn ee ee ee ee 
A. FIELD-EFFECT TRANSISTOR CHARACTERISTIC --------------- 

B. THE TRAC NETWORK ANALYSIS PROGRAM -------------------- 

1. General ----------------..-.-.--.--...-.......-.-.- 

2. Nonlinear Modelling Procedure -------------------- 

III. FIELD-EFFECT TRANSISTOR MODELS --------------------------- 
A. GENERAL -------- occ ec enn eo nn on ne ne ne ee ee ee eee ee 

B. TRIODE REGION ANALYSIS ------------------------------- 

l. Pointwise Linearization -------------------------- 

2. Convergence Method ------------------------------- 

C. PINCH-OFF REGION ANALYSIS ---------------------------- 

l. Pointwise Linearization -------------------------- 

2. Convergence Method ------------------------------- 

D, LUMPED PARAMETERS MODELS ma a a nn ee eee 

1. General ------------------------------~-------+------ 

2. Junction FET Model ------------------------------- 

3. Insulated-gate FET Model ------------------------- 

IV. EXPERIMENTAL RESULTS ------e--- ern ee ne nee nen ene ene eee 
A. PARAMETER DETERMINATION ------------------------------ 

i. Gener | Sea eee eet rarer 

2. Nonlinear Parameters --------------------------+--- 

3. Determination of Capitance Values ---------------- 

B. MODEL PERFORMANCE ------------------------------------ 





ee CONCLUSIONS ---------------- n-ne nr rr nn ee 
A. GENERAL ---------- oor nn nnn nr rr rr ee eee 
B. SUGGESTIONS FOR FURTHER STUDY ------------------------ 
APPENDIX A ------3- rn nn nen nr nn ee en ee ee ee eee 
BLBLIOGRAPHY ----r ncn rnc rrr rrr rr nn ee ee = - 
INITIAL DISTRIBUTION LIST ------------------------- oooncnn------ 


FORM DD 1473 ---------- 2-2-2 ee ee 





I. INTRODUCTION 


When a computer is used to analyze an electronic circuit, the user, 
through the input routine of the analysis program, must provide suffi- 
cient information concerning the circuit parameters and device terminal 
relationships to allow the network equations describing the circuit to 
be determined and solved. This is a simple task for passive linear 
components such as resistors, capacitors, and inductors; but active 
nonlinear devices must be described to the computer analysis program in 
terms of mathematical models which simulate the behavior of the device. 
As field-effect transistors (FETs) have become more common in electronic 
circuits, it has become increasingly important to develop accurate 
models of these devices for use in computer programs used in design and 
analysis of such circuits. 

Small-signal FET models are easily produced for any analysis program 
which includes voltage-controlled current or voltage sources as allow- 
able circuit elements. Any modern basic electronics textbook [1] 
includes a description of such models. These are valid only over a 
small range near the bias point at which the parameters were determined. 
Large-signal models, which are valid over the entire operating range of 
the device, are required for dynamic analysis of logic and switching 
circuits, for very accurate analysis of large-signal "linear" circuits 
(such as various classes of amplifiers, limiters, and mixers), and for 
special circuits in which the FET is used as a voltage-controlled 


current source or as a voltage-variable resistor. 





Today's circuit designer can choose from many computer-aided design 
(CAD) programs which offer a great variety of features including 
frequency response, time response, transfer functions, and pole and 
zero locations; but all general-purpose and most special-purpose analysis 
programs include the basic time-response solution of the network 
equations. ECAP [2] and TRAC [3] are two such programs. ECAP is a 
widely used general-purpose analysis program with time-response output 
only. TRAC is the network analysis program which essed in developing 
the FET models presented in this paper. 

Roberts and Harbourt [4] developed piecewise linear large-signal 
FET models for use with ECAP, ECAP requires that active devices be 
represented by linear passive elements and controlled current sources. 
Branch current-sensing switches are used to change parameter values 
when going from one piecewise linear section to another. ECAP large- 
signal modeling is thus inconvenient, requiring the introduction of 
switches and additional network nodes and elements to generate the 
piecewise approximation to the nonlinear FET characteristics. These 
detract from the total network node and element capability of ECAP. 

The TRAC program is much more flexible and allows the analyst to 
define his own models by writing a FORTRAN subprogram in which he enters 
terms directly into elements of the network matrix equations. This 
paper describes the development of nonlinear models for FETs and the 
procedure for determining the matrix terms using an iterative pointwise 
linearization of the nonlinear FET characteristics. 

The discussion here is given in terms of eee devices although 


FETs of all types are also made with p-type channels which operate in 





the same fashion with current directions and all voltage polarities 
reversed. Inclusion of these complimentary types would have compli- 


cated the discussion without demonstrating anything new or different. 
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II. BACKGROUND 


A. FIELD-EFFECT TRANSISTOR CHARACTERISTICS 

The theory of operation, advantages of using, and techniques of 
manufacturing field-effect transistors (FETs) have been well described 
elsewhere and only an introductory review of their operating character- 
istics is necessary for understanding the computer model of these 
devices. Extensive bibliographies of FET studies [5,6,7] have been 
produced to guide researchers who require more detailed information 
about the devices. 

An FET is basically a piece of lightly doped semiconductor material 
in which the size and shape of the conducting portion, called the 
channel, is modulated by a controlling voltage to prodiice variable 
conductance. The terminal to which the control voltage with respect to 
the channel is applied is called the gate. A salient feature is the 
pinch-off characteristic, in which the channel saturation current is 
limited to a level determined by the controlling gate voltage. 

The channel current is carried by majority carriers under the 
influence of an applied electric field. Mobile carriers travel from 
the end of the channel called the source, to the opposite end, the 
drain, to which the attracting electric field is applied. Many FETs are 
symmetrical, i.e., the source and drain ends of the channel are inter- 
changeable. 

In the junction field-effect transistor (JFET), channel modulation 
is accomplished by depleting its carriers by reverse-biasing of the P-N 


gate-channel junction. This restricts operation of this type to the 





polarity which provides this reverse bias. Since the channel is lightly 
doped and the gate diffusion region is more heavily doped, small reverse 
voltages on the gate are capable of depleting a large portion of the 
channel. 

Insulated-gate field-effect transistors (IGFETs) are surface- 
controlled devices and channel modulation depends on a modification of 
the charge distribution when an electric field is applied to the surface. 
This field is obtained by applying a voltage to an electrode (the gate) 
which is a thin metal film deposited on a solid dielectric material 
insulating it from the channel. The metal is commonly aluminum and the 
insulator is usually an oxide such as silicon dioxide (Si0,). For this 
reason IGFETs are also called MISFETs, for metal-insulator-semiconductor 
FET, and the primary subgroup which uses an oxide insulator is called 
MOSFET, to describe the metal-oxide-semiconductor, structure. 

The channel under the insulated gate is very lightly doped and it is 
possible to apply a field large enough to invert the doping concentra- 
tion. This is used to produce two different types of IGFET. If the 
drain and source contacts are made at regions in which n-type material 
has been diffused, a depletion-enhancement type IGFET, often called 
simply a depletion type, can be made by connecting these two regions 
with a lightly doped n-type channel. By immobilizing charge carriers, a 
negative voltage applied to the gate depletes the channel of usable 
carriers. If the field is high enough, the channel can become intristic 
and inverted to p-type, completely cutting off channel current. If a 
positive voltage is applied to this gate, BIC EOns will be attracted 
into the channel, enhancing its conductivity. Thus this type is capable 


of operating with either polarity of voltage on the controlling gate. 
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Figure 1. Field-effect transistor structures: (a) Junction 
FET, (b) Depletion-type IGFET, (c) Enhancement- 
type IGFET. 
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The third type of FET is the enhancement IGFET. If the drain and 
source n-type regions are diffused into a very lightly p-doped channel, 
conduction will not take place until the channel is inverted by applying 
sufficient positive voltage to the gate. The voltage at which conduction 


begins is called the threshold voltage, V Production of enhancement- 


rT 
type IGFETs requires very light doping of extremely pure semiconductor 
material. Manufacturers are making efforts to reduce the magnitude of 
Ver as much as possible. 

Simplified structure diagrams of these three types of field-effect 
transistors are shown in figure l. 

The IGFET channel is laid out on a lightly doped substrate of 
Opposite doping type. In normal operation the substrate is connected 
to the source, but often the substrate lead is brought out separately 
and a bias voltage or signal may be applied, making the IGFET a four- 
terminal device. Long [8] found that for several transistor types at 
a variety of values of terminal voltage, channel current was nearly 


independent of the substrate bias voltage, V He reported that the 


ss° 
only noticeable effect was (for an enhancement type device) a few 
tenths of a volt variation in Ver for a zero to 20 volt variation in Veg: 
Of interest to the circuit analyst are the terminal I-V relation- 
ships and immitance characteristics of the device. Since the IGFET 
channel current is independent or only weakly dependent on the substrate 
voltage, and since in most applications the substrate and source are 
interconnected, the three types of FETs can be considered as three- 


terminal devices. Dual-gate MOSFETs are true four-terminal devices in 


that channel current is strongly dependent upon all four terminal 
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voltages. Particular device characteristics are dependent upon geome- 
trical structure and placement of the two gate electrodes. 

Early derivations of FET characteristics started by considering 
the impurity profile of the channel and gate-channel contact. This led 
to solutions valid only for the particular doping profile under consid- 
eration. Sevin [9] showed that although Feleeene profiles led to 
expressions of different form, all profiles, including the extreme cases, 
gave approximately the same numerical result. Sevin was thus able to 
choose a profile which gave a simple expression for the square-law 
transfer characteristic. Middlebrook's charge-control approach [10] 
resulted in a simplified derivation valid for all types of FETs and gave 
the same simple expression shown by Sevin. Based on this work, Angelo 
[1] has set forth a set of equations describing the functional dependence 


of channel current on the drain, source, and gate voltages, v and 


yee IS 
Yo respectively. These have been used here. Single-gate FETs have two 
primary operating regions, the triode region in which the device behaves 
as a voltage-variable resistor, and the pinchoff (or pentode) region 
where the current is determined mainly by the gate voltage. The 

boundary between these regions is the locus of points where the gate and 
drain-to-source voltages combine to just pinch off the channel, resulting 


Pa esaturaction. 


In the triode region, 0 < Yong < Yes 7 Vp P 
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(1) 


where Vp is the pinch-off voltage (that voltage which when applied to the 


Patch cealoncectitticient to pinch off the channel) and where Iss is a 
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constant (the drain saturation current) which depends upon channel 


geometry, impurity density and profile, and charge carrier mobility 


(and hence temperature). Inco is defined to be the channel current, ins? 
= d = = P 
at Vos O an ng Vp 
In the pinch-off region, 0<v Se ee oN 


GS R DS ’ 


_ GS .2 
ps pss CV? | 2) 
E 
Here, for an ideal device, the current is independent of Vng? for 
ye eV. - V In an actual device there is a small channel conductance, 


DS GS Ee 


Bap? in the pinch-off region, which varies with v in such a way as to 


GS 
always be small compared to the transconductance effect. 

Equations 1 and 2 are valid for junction FETs and depletion-type 
IGFETs. For 1 < Vac! Vp the channel current is essentially zero. The 
JFET is restricted to operation with the gate biased so as to reverse- 
bias the P-N junction. Although a forward bias destroys the junction, 
transient forward voltages of a few tenths of a volt are permitted and 


equations 1 and 2 apply. 


The enhancement-type MOSFET has a characteristic threshold voltage, 


Vee instead of the pinchoff voltage Vp» but the equations are of the 
same general form: 

jee | |e Vag! Ve : Le 0 (3) 
and for an n-channel device, for 0D<vi._<v - V.., (triode region) 


DS GS ih 


7 Z 

tng = KEW Eg - Vy) Ypg ~ Ypg ? : (4) 
and for 0 < Vas 7 Vip < Vpo (pinch-off region) 

i_. = K(v - V \e (5) 

DS GS Ab ° 
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2 
Here tant ie et i 
K is a constant not equa O Liga! Va since Ings has no meaning in 


the enhancement type device where ine =O whenv,. = 0. K is a function 


GS 


of channel geometry, doping, and effective carrier mobility. 


B. THE TRAC NETWORK ANALYSIS PROGRAM 
1. General 

TRAC (Transient Radiation Analysis by Computer) is a digital 
computer program which simulates the time response of an electronic 
circuit to one or more arbitrary forcing functions. It is a special- 
purpose program in that these forcing functions may be radiation stimuli 
as well as electrical inputs. An extremely flexible program, it is 
suitable for general network analysis in or out of a radiation environ- 
ment, and has features which make it especially suitable for modelling 
of nonlinear devices and components. 

TRAC is capable of initializing and computing the time solution 
of a large general set of simultaneous linear and nonlinear, ordinary 
integral and differential equations, provided the user puts these into a 
finite defference equation of specified form. 

The version of TRAC available for this study was written in 
FORTRAN IV and assembly language adapted for use on the IBM-360/67. 

This version used 158 k bytes of storage. Information given here about 
the TRAC program should be sufficient for understanding the non-standard 
modelling procedure used here. Those who require detailed information 
will find TRAC to be well documented in references 3 and 1l. 

TRAC uses the nodal analysis method which permits topological 
circuit description and data entry. TRAC has built-in models for linear 


circuit elements (resistors, capacitors, and inductors), current and 
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voltage sources, diodes, and bipolar junction transistors (BJTs). It 
also has a routine to automatically write the nodal equations for the 
network using these built-in models with topological input data, and 
solves these equations for the node voltages. TRAC also calculates 
branch currents and powers. 

In TRAC analysis, a network is treated as a set of nodes 
interconnected by network elements as shown in figure 2. The reference 
(ground) node and other independent nodes are not unknown and are not 
included in the analysis. The current ip 1. &t) is defined to be the 

3 
; th , 
current leaving the P node due to the voltage-current (I-V) constraints 
th ; 
on the k network element and the unknown node voltages, at time t. 
This current has a functional relationship to the node voltage, the 


parameters of the network elements, the current and voltage sources, 


and time, t: 


“WO, = eee P pao, (6) 


The exact form of this function depends on the network topology and 
the I-V relationship of the ee element. 
By Kirchhoff's current law (KCL), the total current leaving 


each node is zero, since no charge is stored at the node; therefore 


2 


) = 0 , =‘(for each node P_ . (7) 
k 


ale 


In TRAC, the I-V relationship affecting the current leaving node P for 
the pee element is represented by a linear difference equation which 


can be put into the form 


eo 


ip g(t, =] E Ap g(t) * Yalta) ]- ome gle ® 
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Figure 2. TRAC concept of a set of nodes inter- 
connected by network elements. 
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where P and Q are node numbers, 


o th e 
Ah, Q bt) is the k network element admittance contribution 
> b 


at time t , 
n 
Tia eae is the voltage of node Q at time ce » and 
Ele : 
AW, . 6E) is the k network element current contribution at 
3 
time t_. 
n 
Substituting the current form of equation 8 into equation 7, yields 
another form of Kirchoff's current law: 
) { - bhp g,k’fn) * Yelt,)] - bp (ty) } = : ow 
k 
Writing equation 9 simultaneously for each node P, in the network, 


yields the matrix form: 


H-V = T (10) 
where V is the vector [v_, v v aa! Ae 
le 2? eee Q’ oo NV 9 

See can ; z Ah nvyk So 
H = @ @ . 

i SN ai ae ; : Anny nv, kt? 

B hwe se (€.) = Aw (t_) 

“(i 1,k*on ce: ¢ OVNV, k 0? | 


and NV is the number of dependent nodes. 
The nodal admittance matrix, H, is normally nonsingular and is 
a sparse matrix. Diagonal terms represent element conductances; off- 


diagonal terms occur symmetrically and represent transconductances. 
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TRAC proceedes to solve the matrix equation 10 repeatedly, using 
automatic time step (At) control. An automatic DC initialization is 
only a special case of the transient solution. In the version of TRAC 
available for this study, equation 10 was solved using Gaussian elimi- 
nation with pivoting although the sparseness of the matrix indicates a 
more computationally efficient method could have been used. 

2. Nonlinear Modelling Procedure 

TRAC handles the nonlinearity of the diode primary diffusion 
current equation, which is an exponential function, by pointwise 
linearization and ag modified Newton-Raphson iteration scheme [11]. Using 
this standard diode model, a modified Ebers-Moll transistor model for a 
BJT is constructed using two diode standard models and two controlled 
current sources. TRAC has a capacity for thirty of these diode or 
transistor models to be used in any network. The TRAC User's Guide, 
reference 3, describes how models can be constructed for zener diodes, 
Uiimiunmction Eransistors (UJTs), silicon control rectifiers (SCRs), and 
four-layer PNPN switches. These non-standard models are combinations of 
the TRAC built-in diode, transistor, and linear element models. It is 
significant that no such model is suggested for any type of field-effect 
transistor. 

Any mathematical model of any field-effect transistor must 
involve a square-law transfer function rather than the exponential 
function characteristic of the diffusion current devices. Thus the 
only nonlinear equation built into TRAC is inappropriate for modelling 
of field-effect transistors, and no combination of the TRAC standard 


built-in models will suffice. 
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The TRAC program achieves its remarkable flexibility by calling 
a user-written FORTRAN language subroutine in which the user has access 
to nearly all the program variables during program execution at four 
separate times during each iteration. This is where the user can enter 
any function (such as the square-law transfer function for the FET) and 
the statements for calculating and entering terms directly into the 
network matrix equations. The equations describing the mathematical 
model must be written as Linear difference equations of the form of 


equation 8. The Ah and Aw terms are then computed and entered directly. 
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Til. FIELD-EFFECT TRANSISTOR MODELS 


A. GENERAL 

A good model must posses two important attributes: it must 
accurately simulate the device characteristics over the range of 
intended use, and it must be computationally efficient. It is also 
desireable that required device parameters be economical to obtain, 
either provided on the manufacturer's specification sheet or easily 
measured, 

The pointwise linearization method employed here generates a linear 
small-signal model at each step in the time solution. Of course it is 
not economical to measure device characteristics at each of the possible 
Operating points, nor would storage and retrieval of these parameters 
lend itself to computational efficiency. The parameters appropriate to 
each point are better determined by partial differentiation of the 
device I-V expressions and evaluation of these derivatives at that 
point. If the functions involved are continuous and accurately express 
device behavior, a simple iterative scheme can be devised and conver- 
gence criteria chosen to give a very high degree of accuracy at all 
points. 

Equations 1 and 2 express the I-V relationships for junction FETs 
and depletion-type IGFETs. These equations have the advantage of being 
equal at the locus of points Yong = Yas - Vp ,» which describes the 
boundary between the triode and pinch-off operating regions. The slopes 
of these two equations are both equal to zero at this boundary for an 


ideal device. 
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B. TRIODE REGION ANALYSIS 
1. Pointwise Linearization 


In the triode region, 0 < Yns < Vag Vp ,» partial differenti- 


ation of equation l gives 


i I 
DS DSS | | 
= =- 2 pS i 


At an arbitrary point (VDS, IDS) on the output characteristics curves, 





figure 3, the pointwise linear approximation to the output curve is the 


straight line 


ieeeeus (ve. - V0s) + IDS (12) 


where GDS is the slope at (VDS, IDS) given by 


; _ {pss ; : 
GDS = ae ar acves eee ws | (13) 


(VDS, VGS) uD 


Use of FORTRAN notation for TRAC variables from this point is 


more convenient and adds clarity. In TRAC notation, the array terms are 


denoted by 


5 ye ae H(P, Q) and 
AW, 468) ek) 


Putting equation 12 into the form of equation 8 and using FORTRAN 


notation for the variables, gives 


ing = (GDS) v, (t_) + (-GDS) v, (t,) - (GDS * VDS - IDS) 
Since es one = Toe (S and poo sne - si vatc a) 


a 





boundary 


ee 


“ps - P 
V 


PiIrseetcial solution 


i 
DS = 0 
GS 
point (Vag Ys ) / 
I 1 / 
slope = GDS 
Pi Yes — VGS, 
IDS, = —=— —— ; / 
slope = GD 
| / : 
| / 
i 
p ) / 
0 | VA 
IDS, a8 : J 
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| /f 
4 
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oe load line 


V 
DS 
(O70) VDS VDS 


Figure 3. Triode region convergence procedure showing 
One iteration from initial point Py to the 
PeGtmtartal solution point, Pi. 
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-GDS 


H(D, S) 
T(D) = GDS * VDS - IDS 


-GDS 


and HCS=aD) 
H(S, S) = GDS 
T(S) = IDS - GDS * VDS 
where D and S are the network node numbers to which the drain and 
source, respectively, of the io field-effect transistor are connected. 
2. Convergence Method 
Choice of the starting point, (VDSp» IDS) Mt Leure 3 eis 
arbitrary but the TRAC DC initialization program begins with all node 
voltages set equal to zero. During the time solution it is both 
logical and convenient to start each new time step C. using the voltages 


v(t D calculated at the previous time, t at which a solution was 


n-1’ 
accepted. These voltages are automatically updated by TRAC each time 
a new time solution is accepted and time is advanced by the current 


value of At. 


For each iteration within a time step, VGS and VDS are given 


by 
= - 14 
MOS 2 (14) 
and VDS = Vn ~ Ve (15) 
where Vor Yyp and Vo are the current, or trial solution values, of the 


FET terminal voltages. GDS is given by equation 13. Then let 


iE 
IDS 7> 228 2(VGS - V) VDSe5 (vps)" | (16) 
V 
le 


25 





Using these values, the matrix terms described above are then entered 
and the network matrix equation, equation 10, is solved for a trial 


solution, 
Vee eee Vi, OV, Vy oy V Ar (17) 
i este ete oee GC” > “NV 


The load line shown in figure 3 represents the linearization of 
the I-V relationship of the equivalent circuit of the external network 
looking in at the nodes to which the source and drain are connected. 
The only thing known about this line is that if the point (VDS, IDS) 
satisfies the network matrix equation, then this load line will pass 
through (VDS, IDS). In any case, this load line will pass through the 
feeetetrial solution point, (VDS, » IDS,); which is determined by 
substituting the trial solution into equations 15 and 16. This 
represents a return to the output curve for v_. = VGS, along a line 


GS 1 


of constant voltage, VDS, This new point is used to calculate 


(VDS IDS,),; the next iteration, by the same procedure. This process 


Me 


is iterated until the convergence criteria, 


| vDs - Vy =) vr <S § and 
cs vt ve| = 6 


are satisfied, where § is a small, arbitrary constant. At this time the 
trial solution is accepted, and the TRAC program updates the program 


variables and time is incremented by At to th 


41° The above procedure 


is then repeated for this next time step. 
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C. PINCH-OFF REGION ANALYSIS 
1. Pointwise Linearization 


In the pinch-off region, 0 <v ee eV for an ideal 


Gs -p~ ‘ps’ “ps 


device is given by equation 2. In an actual FET, Lic in the pinch-off 


region is not completely independent of v The magnitude of this 


DS’ 


effect depends on the design of the transistor. The gradual increase 


in Ins with increasing Vive is partly due to the fact that channel 


length does not remain constant under pinch-off conditions as was 
assumed in the theoretical analysis leading to equation 2, but decreases 
slightly with increased voltage. This has a greater effect in short- 
channel devices designed for high-frequency applications. [It can be 
neglected in audio FETs. This non-zero slope, By can easily be mea- 
sured at any bias point, and for a small-signal model, incremental 


drain current is given by 


: ze i 
*ds En Yes bq Yds 


For large-signal modelling by pointwise linearization, the functional 


relationship i f(v ), must be known in analytic form. The 


Da Ese DS 


partial derivative di) 6/O%p> can then be evaluated exactly as it was in 
the triode region case. This functional relationship is not known in 
general but can be approximated empirically for a given device by mea- 
surements, if the desired degree of accuracy so requires. 


Since in equation 2 i is independent of V 


DS? Ens is zero, and 


DS 


we let 


2 = S| ee cs 
veces "p op 





(18) 
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Figure 4. 





Pinch-off region transfer curve showing 
pointwise linearization procedure. 
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This is the slope of the transfer characteristic curve, figure 4. At a 


point (VGS, IDS) the tangent line is given by 


ecg) VGS) IDS (19) 


where 2 
on Dss [1 : aa (20) 


Ny 
V 


GS = VGS 
Putting equation 19 into the form of equation 8, 


Lge? = (GM) Ylen) + (-GM) CRY - (GM * VGS - IDS) . 


As in the triode region, Lye) = ing ft) and eee an = - ing(t) 
Therefore 
H(D, G) = GM 


H(D, S) = -GM 
H(S, G) = -CM 
H(S, S) = GM 
T(D) = GM * VGS - IDS 


~GM * VGS + IDS 


T(S) 
where D, S, and G are the network node numbers to which the drain, 


: th 
source, and gate, respectively, of the k FET, are connected. 


2. Convergence Method 
Convergence is achieved in the pinch-off region by the same 
method employed in the triode region. For each iteration within the 


Gine Step, t. 


VGS = Ya Yo 5 C21) 


and GM is calculated by evaluation of equation 20. Then IDS is given by 


b 
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IDS = Lee [ - ves P ; (22) 


Starting at this point, (VGS, IDS), the above values are used to 
determine the matrix terms in equation 10 which is solved for a trial 
solution, 


iE 


anens Vir : 


V= [vy> Vos ses Vos Vos Yoo 


The trial solution value Vos is compared with the value of VGS to 
determine if convergence has been achieved, in which case the 

solution is accepted and program variables updated and time incremented. 
If convergence is not achieved, the solution proceeds back to equation 


21 using the trial solution values to determine the next point 


(VGS, ; IDS, ) and a new value for GM. 


D. LUMPED-PARAMETER MODELS 

1. General 

A physical device such as an FET operates by propagation of 

fields and movement of charge carriers between points in space. This 
propagation and motion takes time which results in delays in trans- 
mission of information between the device terminals. The interaction 
of fields and charge carriers in an FET is distributed throughout the 
entire active channel under the gate. Thus exact simulation of the 
behavior of this device requires a model in which both time and space 
variables are included. The models described in this section do not 
provide time delay and lumped analysis is used. 


Reference 12 describes a distributed IGFET model and an 


associated computer program which includes space variables as well as 


28 





time to simulate IGFET behavior to high precision. This model is valu- 
able for analysis of extrinsic effects and as an aid to design by 
prediction of the characteristics of proposed FET designs. This type 
of analysis model is incompatable with available second-generation 
network analysis programs since these programs require that distributed 
effects be lumped into discrete network elements. 

In figure 1 it can be seen that the active channel under the 
gate in each of the three types of FET has a voltage drop distributed 
along its entire length. Thus voltage-dependent capacitance and 
leakage effects are distributed and continuously variable along the 
gate-to-channel interface. 

2. Junction FET Model 

Consider how the JFET gate appears to an external circuit. The 
channel depletion effect is distributed along the length of the gate 
P-N junction. This junction is reverse-biased and appears as a 
distributed capacitance with a high shunt resistance. The device 
control voltage is that which appears across this capacitance and is not 
necessarily equal to the input voltage due to the time required to 
discharge or charge this capacitance through the gate bulk semiconductor 
resistance. This capacitance must be split between the drain end of 
the channel and the source end to get the required lumped model. It is 
not equally distributed since the drain and the source are biased 
differently with respect to the gate in normal operation. The portion 
of the junction near the drain end of the channel is normally strongly 


2 


reverse-biased and the capacitance value here is nearly independent of 
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Figure 5. Lumped field-effect transistor models. 
(a) Junction FET, (b) Insulated-gate 
FEL, 
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Vv This capacitance and the associated shunt conductance can be 


GS’ 
modelled using the TRAC standard capacitor model which includes pro- 
vision for a shunt resistor. 

The source end of the junction is less strongly reverse-biased 
or even slightly forward-biased, depending on the instantaneous value 


of v It exhibits the voltage-variable capacitance and conductance 


GS’ 
properties of an ordinary P-N junction diode. A lumped model of this 
gate-to-source junction can be produced by using the TRAC standard 
built-in model for a junction diode to simulate these effects. 

These lumped approximations are shown in figure 5(a). The 
drain-source circuit has no current paths to the gate except for the 
leakage and capacitve effects mentioned above. Transconductance and 
output conductance are simulated by the nonlinear current source, Lis: 
which is modelled by pointwise linearization as described earlier. The 


drain circuit capacitance is lumped into G and is represented by the 


ds 


TRAC standard capacitor model. 
3. Insulated-Gate FET Lumped Model 
Figure 1 shows that the gate of an IGFET resembles a parallel- 
plate capacitor with a thin dielectric layer. The dielectric results 
in a very low gate conductance, and since the thickness is constant, 
the distributed capacitance is uniform along the channel and independent 


of v This capacitance can be split between the source and drain 


GS° 
ends of the channel and both lumped capacitors simulated by the TRAC 
standard capacitor model. This approximation is shown in figure 5(b). 
This model can be used for both depletion-type and enhancement-type 


IGFETs, the only difference being in the expressions used to obtain the 


linearization of Ins? Equations 4 and 5 which are applicable to the 
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enhancement-type IGFET are of the same form as equations 1 and 2 
respectively, for the depletion-type IGFET. Thus the pointwise lineari- 
zation and iteration process for convergence is the same for both types 


ot IGFET, 
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IV. EXPERIMENTAL RESULTS 


A. PARAMETER DETERMINATION 
1. General 

Accurate simulation of FET circuits requires not only a model 
of suitable form but also accurate parameter data to be used with the 
model. It would be ideal if the required information could be obtained 
from the manufacturer's specification sheet for a particular device. 
With present fabrication techniques useable FETs of a particular type 
are produced with a wide spread of parameter values; to demand close 
tolerances would result in low yields. Therefore specification sheets 
list a range or "typical" parameter value, and where a particular 
parameter is critical to some aspect of performance, a maximum or mini- 
mum value is listed. This insures that any given device of that type 
will meet certain minimum performance criteria although many will exceed 
them. 

Often the circuit being analysed must be insensitive to para- 
meter variations and must perform satisfactorily under worse-case 
conditions. In such a case, the maximum, minimum and typical values 
given on the specification sheet can be used as a starting point for 
worst-case and sensitivity analysis. In other cases the circuit to be 
simulated will be one in which the circuit is to be carefully adjusted 
for optimum performance with a particular device installed. In this 
case all critical parameters must be carefully measured to insure 
accurate simulation. Reference 4 describes methods of measuring FET 
parameters. Noncritical values may be taken from the specification 


sheet, if available, or estimated. 
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2. Nonlinear Element Parameters 
The pinch-off voltage, Vp and the drain saturation current, 


IL are two of the parameters required in the models presented here. 


DSS’ 


While Ings is very easily measured it is difficult to measure Vp with 


much accuracy, since the transconductance and current become very small 
as Vp is approached. We can determine Vp much more accurately by mea- 


suring a different parameter. For the pinch-off region, the transcon- 


ductance is given by equation 18, 


_ *Iygs Met: 
es ey 
P Ie 


By evaluating g. at v 5 = O where g is maximum, 


G 


_ *Ihgg 





8 aaa ; 
max Vp 
Then ZY 
V - . —D55 y (23) 
P g 
max 
Since gnax is easily determined at the same bias point as Ings they can 


be measured simultaneously. The circuit used in this study to measure 

these values is shown in figure 6. Table I lists parameters determined 

by this method for two types of FETs chosen for this study. Accuracy is 
TABLE I 


Parameter values for two types of field-effect transistors 


Hee hss Emax Vp ee eae as 


Type (mA) (umho) (volts) (pf) (pf) (pf) 


2N3819 (JFET) 4.1 4000 -2.05 8. 4, A 
40468A (IGFET) 13.2 6010 ~-4.39 5.38 0.12 1.28 
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Figure 7. 


Pulse inverter circuit. Node numbers 
used in the TRAC program are encircled. 
EI(k) is the FORTRAN coding for the k* 
grounded voltage source. The friodes G, 
Sl, and S2 are not independent nodes and 
do not enter into the set of simultaneous 
equations solved by TRAC analysis. 
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indicated by the number of significant figures used. Results must be 
interpreted with caution since the computer program expresses numerical 
results to many more significant figures than can be justified on the 
basis of the model and accuracy of the empirically determined parameters. 
In modelling the junction FET, a TRAC built-in diode model has 


been used. This model requires that the parameter I Ss be supplied. 


G 


Tags may be determined by shorting the drain to the source, applying a 


reverse bias, and measuring the source-to-gate current under this 
condition. If the FET is never to be operated in a forward-biased 
condition, this parameter is not critical since the current will be 


insignificant, and the value may be taken from the specification sheet. 


3. Determination of Capacitance Values 


The models presented in figure 5 have three capacitance values 


The TRAC diode model used 


which must be determined: C ,C, , and C, . 
gs ds 


dg 


the maximum at v = 0 of 


in the JFET also requires a value for C30? Gs 


the voltage-variable diode capacitance. This is the bias level at 
which Ce ss and Cae are measured. The input capacitance, Cs ss? the 
output capacitance, Crane and the reverse-transfer capacitance, Cae? 
were very small for the two FETs used in this study. While it would 
have been possible to measure these values, they were taken from the 
specification sheets. The capacitances required in the model can be 
determined from these parameters. Since C. os is the input capacitance 
with the drain shorted to the source, C. os = oe + oie Similarily, 


C Cet C and C = aC 
dg 18 


» These relationships can be solved for 
oss ds SS dg 


the required values: 
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. 


dg rss 

C = C, = 1S 

gs iss rss 
& = C - C 

ds oss rss 


Values determined for these capacitors are shown in table I. In the 
FET, ve is supplied to the diode model, which automatically varies 


gate capacitance with voltage. 


B. MODEL PERFORMANCE 

Large-signal dynamic models, while capable of simulating the opera- 
tion of small-signal amplifiers and other small-signal circuits, are 
more importantly capable of simulating switching transients. For this 
reason the simple pulse-inverter circuit shown in figure 7 was chosen 
to demonstrate the ability of the proposed models to simulate device 
behavior accurately. This circuit was coded for TRAC analysis and was 
physically constructed with the same FETs which had been used in para- 
meter determinations. Waveforms were observed and signal levels 
recorded to provide a basis for evaluation of model performance. 

Figure 8 shows the junction FET model response to a 0.35-microsecond 
pulse with a rise time of 5 picoseconds. The input pulse level went from 
V.. to zero to turn on the device. Comparison of model transient response 
to actual transient response is summarized in table II. 

The acutal device was pulsed by a Datapulse type-101 pulse 
generator which has a pay herr eee rise time. The voltages v and 


GS 


Yng were displayed on a Tektronix type-551 dual trace oscilloscope 


triggered by the pulse generator. A type 53/54K preamp having a rise 


Oy) 
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Figure 8. Junction FET model response to a 0.35 micro- 
second input pulse of a 5 nanosecond rise time. 


= 





time of 0.006 microseconds was installed. Measurements were made with 
compensated RF probes but no other attempt was made to reduce stray 
Capacitance in the circuit. 

The computer programming for the model of this circuit is included 
in Appendix A. Subroutine TRAEQ is called by TRAC to supply the re- 


quired matrix terms for the nonlinear portion of the model. 


TABLE II 


Comparison of switching levels and transients 
for FET model and actual device 


ON Rise OFF Fall 
delay Time delay Time 
(volts) (volts) CSCC) MC ceGy (see) mm (sec) 


PEL Vg (ON) Vig (OFF) 
Ty pe 


2N3819 JFET 


actual orice: Loyal O20TS 0.102 3 0n1 0.100 
model Ste) 14.99 0.0135. 502 055an 0. 0U> 0.075 
40468A IGFET 

actual Liste! 15.0 OTO16 S050 20 mee0 2005 0.070 
model 1.40 14.79 Oc 0500671050 0.005 





Figure 9 shows the IGFET model response to a similar exitation. 
The input pulse level was changed to match Vp for the IGFET. Perfor- 
mance comparison with the actual circuit is included in table II. Rise 
time was much faster with this FET since it is a high-frequency type. 
Measurements were made with the same equipment described above. Sub- 
routine TRAEQ was the same as that used for the JFET. The TRAC data 
deck was modified to enter data Appropriate to the type 40468A 
parameters. 

Table II shows that in most cases the compared values were reason- 


ably well in agreement. The ON voltage of the IGFET was substantially 
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lower in the model than in the actual device. This was attributed to 
the determination of Ino and ae by using a CW signal at zero gate 
bias which generated quite a lot of heat. This resulted in using a 
value for Ing in the model which was valid for a temperature higher 
than was actually present during the low-duty-cycle pulsed operation of 
the inverter circuit. Even so, the value pened was well within the 
expected parameter spread for this device type. 

Actual rise and fall times were consistantly greater than simulated 
values. This was to be expected since circuit stray capacitance was 
not included in the model circuit. Except when the FET is part of an 
integrated circuit, stray capacitance will be significant compared to 
device internal capacitances, and should be included where rise and fall 
times are critical. 

The delay times in the JFET model response were due to the time 
constant associated with the charge storage capacitance which is in- 
cluded in the TRAC diode model. This resulted in a short spike which 
can be seen in figure 8 at the switch-on point. Actual circuit delay 
times were quite small. Delay times in the IGFET model response were 
zero since no delay was included in the model. 

TRAC permits as many as 200 iterations before signaling nonconver- 
gence of nonlinear models. Observation of intermediate results during 
the TRAC solution revealed that usually not more than two iterations 


were required to achieve convergence and quite often the initial solu- 


tion was accepted. 
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V. CONCLUSIONS 


A. GENERAL 

The field-effect transistor models proposed here simulated actual 
device behavior quite well over the normal large-signal operating range 
associated with switching circuits. For a reasonably accurate simula- 
tion the parameters required for the model were few and easily deter- 
mined. The models are fully compatable with TRAC, a very flexible 
network analysis program, and converged rapidly, as was expected, result- 
ing in computational efficiency. These models lack a time-delay feature 
but otherwise have excellent transient behavior. 

Emphasis was placed on empirical determination of parameters; how- 
ever these might also be determined from the material characteristics 


and structural geometry of the device being modelled. 


B. SUGGESTIONS FOR FURTHER STUDY 
1. Simulation of Multiple-Device Circuits 

While the circuit chosen to demonstrate the models was a single- 
device network, the matrix terms which were derived are applicable to 
any number of devices in the same network. The expressions were given 
in terms of the general gate node, source node, and drain node. Where 
more than one FET is included in a circuit the appropriate network node 
numbers must be assigned for each terminal. Where leads from more than 
one FET are connected to the same set of nodes the more general expres- 
sion for adding the mas DE Leconber olUttoneic ence Oy — HCP QO) + HCP. VE: 
No coding problems should arise but convergence may be a problem in 
circuits with feedback and interaction between several device terminal 


voltages. 
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2. Dual-gate IGFET Simulation 


Drain current in a dual-gate IGFET is a function of Ys and the 
voltages on the two gates. An approximate analytical expression for 
this functional relationship must be obtained to implement pointwise 
linearization. One of the problems anticipated is that there is no sim- 
ple division into operating regions as in the case of the single-gate 
IGFET. Since there are two independent gate electrodes there are several 
operating conditions possible at each point along the channel; further- 
more, these conditions can exist simultaneously at different points in 
the same structure. A particular set of boundary conditions applies to 
the continuity equation for each possible combination of channel condi- 
tions. 


3. Temperature Effects 


a V_, and Ve are related to carrier mobility which is tem- 


DSS’ P? 
perature dependent. TRAC includes temperature effects in its diode and 
bipolar junction transistor models and the temperature can be varied as 
a function of time and network variables by statements in the TRAEQ 
subroutine deck. Given expressions for parameter temperature dependence 
it should be simple to extend these models to include temperature 
Ci reciuc.. 

4, Other Factors Affecting Model Accuracy 

The models presented here assumed negligible output conductance 

in the pinch-off region. Further work is required to develop an approx- 
imate analytical expression for this conductance which can be linearized 
to provide a more accurate solution in the pinch-off region. 


If a FET is to be operated outside of its normal range, it is 


necessary to include breakdown effects. It is possible that these 
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effects could be modelled using standard TRAC elements. The possibility 
of handling destructive breakdown by logical programming in the TRAEQ 
subroutine can be investigated. 

In an actual device the active channel must be extended beyond 
the gate area to permit the source and drain contacts to be attached. 
Since the channel is lightly doped this adds significant bulk resistance 
to each of the terminals of the device. Further, there are small 
contact potentials which have been neglected. These and other extrinsic 


effects must be included for precision simulation. 
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MODELS 


TRAC CALLS THIS SUBROUTINE TO GET MATRIX TERMS FOR NONLINEAR 
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THE FOLLOWING CARDS CONSTITUTE THE DATA DECK FOR THE 2N3819 JFET MODEL 
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